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The conformal Einstein equations and the representation of spatial infinity as a cylinder 
introduced by Friedrich are used to analyse the behaviour of the gravitational field near 
null and spatial infinity for the development of data which are asymptotically Euclidean, 
conformally flat and time asymmetric. Our analysis allows for initial data whose second 
fundamental form is more general than the one given by the standard Bowen-York Ansatz. 
The Conformal Einstein equations imply, upon evaluation on the cylinder at spatial infinity, 
a hierarchy of transport equations which can be used to calculate asymptotic expansions for 

O" 1 the gravitational field in a recursive way. It is found that the solutions to these transport 

equations develop logarithmic divergences at certain critical sets where null infinity meets 

(3JT), spatial infinity. Associated to these, there is a series of quantities expressible in terms of the 

initial data (obstructions), which if zero, preclude the appearance of some of the logarithmic 
. £^ ' divergences. The obstructions are, in general, time asymmetric. That is, the obstructions at 

the intersection of future null infinity with spatial infinity are in general different from those 

5_l ' obtained at the intersection of past null infinity with spatial infinity. The latter allows for 

the possibility of having spacetimes where future and past null infinity have different degrees 
of smoothness. Finally, it is shown that if both sets of obstructions vanish up to a certain 
order, then the initial data has to be asymptotically Schwarzschildean to some degree. 

Pacs: 04.20.Ha, 04.20.Ex 



1 Introduction 

The regular initial value problem near spatial infinity introduced by Friedrich in 15 has proved 
an invaluable tool for the understanding of the behaviour of the gravitational field in the region 
of spacetime near spatial and null infinities. This setting, relying entirely on the properties of 
conformal structure, is such that the equations and data are regular with spatial and null infinities 
having a finite representation with their structure and location known a priori. In order to keep 
the amount of computations at bay, Friedrich's original analysis was restricted to time symmetric 
initial data sets which admit a smooth conformal compactification at infinity. Central to his 
discussion, laid a representation of spatial infinity as a cylinder on which the Conformal Einstein 
equations fully reduce to a set of transport equations. The crucial feature of these equations is 
that through its solutions it is possible to relate properties of the spacetime at null infinity with 
properties of the initial data sets, effectively allowing us to identify the parts of the initial data 
responsible for the non-smoothness of null infinity. In particular, in |15j a regularity condition 
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was formulated, that if satisfied, would preclude the appearance of a certain type of logarithmic 
divergences in the solutions of the transport equations at the sets where null infinity intersects 
spatial infinity 1 . Because of the underlying hyperbolic nature of the Conformal Einstein equations, 
it is to be expected that the logarithmic divergences would propagate along the generators of null 
infinity and leave an imprint on various radiative properties of the spacetime, the decay of the 
components of the Weyl tensor among them. 

The solutions of the transport equations at spatial infinity can be calculated by means of what 
is essentially a decomposition of the unknowns in terms of spherical harmonics. The feasibility 
of this approach was explored in |19| . Leaving aside the calculational complexities, this approach 
is completely algorithmic, and hence, entirely amenable to an implementation in a computer 
algebra system. In {24] this idea was put in practice for conformally flat initial data sets which 
are conformally flat. This family of initial data sets would satisfy trivially the regularity condition 
found in The calculations performed by means of scripts written in the computer algebra 
system Maple V revealed that, generically, further logarithmic divergences would arise. Associated 
with these divergences there is a hierarchy of quantities written entirely in terms of the initial data 
(obstructions), that, if set to zero, would eliminate some of the logarithmic divergences from the 
solutions of the transport equations. Furthermore, an initial data set for which the obstructions 
vanish up to a certain order, would be found to be asymptotic Schwarzschildean to a certain 
order s* which would depend on p» — i.e. s* = s*(p*). This leads to conjecture that the only 
time symmetric, conformally flat initial data sets yielding a development admitting a conformal 
compactification are the Schwarzschildean ones. 

Stationary (vacuum) spacetimes, and in particular the static ones, are known to admit a 
smooth null infinity — see ^J 2 - ^ n order to be able to consider static solutions one has to move 
away from the class of conformally flat initial data sets as the Schwarzschild solution is the only 
static solution having slices which are conformally flat near infinity. In |2.'S| the calculations 
carried out in were extended to the case of non-conformally flat spacetimes which satisfy 
Friedrich's regularity condition. As expected, the calculations rendered a generalised hierarchy 
of obstructions, which, if satisfied, imply asymptotic stationarity to a certain order. Thus, it 
seems that time symmetric initial data sets would render developments admitting a null infinity 
if and only if they are asymptotic static. Providing a proof of this conjecture would constitute a 
remarkable feat. 

Why so much interest in spacetimes admitting a smooth null infinity? The driving force of the 
programme initiated in |15j is to verify the so-called Penrose proposal. Loosely speaking, Penrose's 
proposal states that the gravitational field of spacetimes describing isolated systems — e.g. two 
coalescing black holes — should satisfy in its asymptotic region the Peeling behaviour. That is, 
the components ^ n of the Weyl tensor (where here we are using the standard Newman-Penrose 
notation) should decay along the generators of future oriented light cones as 

*n = o(^r), n = 0,...4, (1) 

where r is an affine parameter along the generators of the light cones. As it was already pointed 
out in , spacetimes with obstructions like the ones found in |^] [53] cannot be expected to 
peel. Thus, the question is now how restrictive is Penrose's proposal when discussing systems of 
physical interest? In any case, it should be noted that non-static spacetimes satisfying the peeling 
behaviour do exist. These have been constructed by Chrusciel & Delay |S] from initial data sets 
which are asymptotic Schwarzschildean outside a compact set. The initial data sets have been 
obtained by using a refinement of the deformation techniques firstly introduced by Corvino in |H] . 
Perhaps one day numerical simulations will allow us to contrast the physical differences (if any) 
of a spacetime evolved, say, from Bowcn-York initial data — which shall not peel — and those of 

lr This condition is expressed in terms of the Cotton- York (Bach) tensor and its derivatives to all orders. More 
precisely, 

S(D n ■ ■ ■ D lq B jk ){{) = for g = 0, 1, . . . , 

where D; denotes the covariant derivative of the metric on the initial hypersurface, Bj k is its Cotton tensor, and 
<S denotes the operation of taking the symmetric, tracefree part of a tensor. 
2 Note, however, that the original proof had a gap which was filled in 7 . 
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the time development of a data and initial data set which exactly is Bowen-York inside a compact 
set, and exactly Kerrian in the asymptotic region — which, arguably, should peel. 

The present article and a companion one — which we shall refer as to paper II — extend the 
analysis initiated in 1241 l2"3"] to the case of non-time symmetric (or time asymmetric) initial 
data sets. In particular, the present article is concerned with those time asymmetric initial data 
sets which are conformally flat. Thus, it encompasses, but it is not limited to, some of the initial 
data sets which are routinely used in the numeric simulations of spacctimes with black holes like 
the Bowen-York and Brand-Brugmann We shall consider second fundamental forms with 

higher order multipole moments like the one used in [Hi]- On the other hand, paper II shall 
deal with time asymmetric initial data sets which are non-conformally flat, like the one that can 
be derived from the Kerr metric written in Boyer-Lindquist coordinates. Most of the ideas and 
methods to be employed in paper I will be developed and discussed in the present article. 

A computer algebra based approach, similar to the one used in |24II23| will reveal the existence 
of a hierarchy of obstructions to the smoothness of null infinity. In contrast to the time symmetric 
case where the vanishing of the obstructions implied a certain degree of regularity at both future 
and past null infinity, the obstructions arising when contemplating non-time symmetric initial 
data can be time asymmetric. That is, there is a hierarchy of obstructions ensuring the smoothness 
of and another one for J 5 " ~ , and these do not imply each other in general. This novel feature 
opens the possibility of the existence of spacetimes where the two disconnected parts of null 
infinity have different degrees of regularity. The existence of this kind of situations had been 
speculated long time ago by Walker & Will |2"fjl |2~7] 3 . 

The analysis undertaken in |201 1251 l2~3"| has shown that the only stationary spacetime that 
can be expected to admit conformally flat slices is the Schwarzschild spacetime. Thus, based on 
prior experience, one would expect that the new hierarchy of obstructions that are obtained in 
this analysis should, if vanishing, restrict the initial data set to be asymptotic Schwarzschildean 
to a certain order. This is found to be the case only if both the obstructions at ^ + and J?~ 
vanish. The Schwarzschild data referred to in this context is not the standard time symmetric 
t = constant slice obtained from writing the Schwarzschild metric in isotropic coordinates, but a 
non-time symmetric, conformally flat slice having a second fundamental form decaying like 



as \y\ — > oo in some coordinates adapted to the asymptotic region — see below. Whether it is 
possible to have a spacetime whose future null infinity is fully smooth, without having to restrict 
the behaviour of the field at past null infinity is a question which would require a substantial 
refinement of the computer algebra methods used in this article, and which will be left for the 
future. 

This article follows, in as much as it possible, the notation and nomenclature of |15| . and in 
particular that of |18| . It makes use of a number of results and techniques which although they 
are available in the literature, they may not be, widely known. Therefore, some time is spend in 
introducing the necessary ideas and concepts. The reader is in any case remitted to the original 
articles for full details. 

The article is structured as follows: in section 2, a discussion of the constraint equations in 
the light of our intended applications is performed. In particular, some results due to Dain & 
Friedrich, guaranteeing the existence of the type of initial data under consideration are recalled. 
In section 3, a discussion of the description of the geometry of the initial hypcrsurface near 
spatial infinity in terms of spatial spinors is given. In section 4, a discussion of the solutions to 
the conformally flat momentum constraint within the space spinor formalism is given. Section 5 
does likewise with the Hamiltonian constraint. Section 6 is concerned with the particularities of 
the propagation equations implied by the Conformal Einstein equations in the region of spacetime 

3 In these two articles, the authors discussed by means of post-Newtonian methods the effects of past null 
infinity of the incoming gravitational radiation produced by two gravitating bodies coming from infinity, i~, on 
approximately hyperbolic orbits which then escape to infinity, i+. 
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near null and spatial infinity. Section 7 deals with some properties of the asymptotic expansions 
near null and spatial infinity introduced by Friedrich. Section 8 discusses the results of our 
computer algebra calculations. Finally, in section 9 some concluding remarks are given. 

In order to perform the calculations here described, a number of assumptions have been made. 
In order to help the reader to keep track of them, they have been clearly marked in the text and 
numbered from 1 to 5. 



1.1 General conventions 

In this work we shall be considering spacetimes (.M, arising as the development of some 
Cauchy initial data (5, h a p,Xafi)- Tilded quantities will refer to quantities in the physical space- 
time, whereas untilded ones will denote generically quantities on an unphysical — i.e. conformally 
rescaled — spacetime. The indices /i, A, . . . (second half of the Greek alphabet) are spacetime 
indices taking the values 0, . . . , 3; while a, /3, 7, . . . are spatial ones with range 1, . . . , 3. The latin 
indices a, b, c, . . . will be used in spatial expressions which are valid for a particular coordinate 
system (usually a Cartesian normal one) and take the values 1, . . . , 3. The indices, i, j, fc, . . . 
are spatial frame indices ranging 1, . . . , 3, while A, B, C, . . . will be spinorial indices taking the 
values 0, 1. Because of the use of spinors, the signature of g^ v will be taken to be (+, — , — , — ), 
and the 3-dimensional metric h a p will be negative definite. Under our conventions and labeling, 
the Einstein constraint equations read 

r-x 2 + Xo/?X Q/3 = 0, (3a) 
D a x a p - D P x = 0, (3b) 

where D denotes the connection associated with the 3-metric h a p, r is the corresponding Ricci 
scalar and we have written \ = X a a - 

The initial data set (S, h a p, Xafi) is to be assumed asymptotically flat in the sense that there 
exists a compact subset of S such that its complement is the union of disjoint subsets Sk, k — 
1, 2, . . . , K of <S — the asymptotically flat ends — each of which can be identified diffcomorphically 
with {R 3 I \y\ > ro} where ro is some positive real number, and |y| = \J (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 3 - 
We shall require that the coordinates so introduced satisfy 



l-ab 



l + ^W + OlrV), as \v\ ~~ * 00 ! m k = constant. 

\y\ J \\y\ 2 J 



(4) 



Moreover, we shall require the second fundamental form to decay as 



* ab = ° ( |^p) ' aS l y ' ^ °°' ( 5 ) 

Besides the above asymptotic flatness requirements, we shall further assume (S, hap) to be asymp- 
totically Euclidean and regular in the sense that there exists a 3-dimensional, orientable, smooth, 
compact Riemannian manifold (5, h a (j) with points S <S, k = 1, 2, . . . q, q some positive integer, 
a diffeomorphism $ of S \ . . . , i q } onto S and a function ft £ C 2 (S) n C co (S \ . . . , i q }) 
with the following properties: 

O = 0, dd — 0, Hess(ri) Q/ 3 = — 2/i Q/ 3, (6a) 
fi>0on«S\{ii,...,i 9 }, (6b) 
h afi = Sl 2 $»/i Q(3 on S \ {ii, . . . , i q }. (6c) 

Note that suitable punctured neighbourhoods of the points ik correspond to the asymptotically 
flat ends of (S, h a p). Thus, each point ik represents a spatial infinity. The discussion in this article 
will be concerned with the behaviour of the gravitational field in one of these spatial infinities. 
So, without loss of generality we assume that there is one of them, which we denote by i. The 
manifold S will be identified via $ with <S \ {i}. 
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2 The conformal constraint equations 

In the sequel, we shall consider asymptotic expansions near infinity for a wide class of solutions 
of the constraint equations Ij^ajl and l|3bjl . Instead of working in the physical spacetime manifold, 
we shall consider a conformally rescaled version thereof 

g^u = ^ 2 g^v (7) 

The first and second fundamental forms induced by the metrics g^ v and g^ on S are then related 
via 

h a p — tt 2 h al 3, Xaf3 = 0(Xa/3 +T,h a ff), (8) 

where E denotes the derivative of Q along the future directed g-unit normal of S. As a consequence 
of the latter two equations the traces x = h a,3 Xaf3, X — h al3 x.a/3 satisfy the relation 

n x = X + 3E. (9) 

The constraint equations (|3a[) and (|3 h>|> can be written in terms of the conformal fields as — see 

na- 

2VLD a D a tt - 3D a ilD a n + ^n 2 r - 3£ 2 - ift 2 - - ( X 2 - X^X^) + 2fffix = 0, (10a) 
n 3 D a (n- 2 Xa p) - V (D f jx - 2VL- 1 D (j Y) = 0, (10b) 

where D a denotes the Levi-Civita connection and r the Ricci scalar of the metric h a /3- The latter 
expressions lead to our first 

Assumption 1 (maximal initial data). In this work we shall assume that 

E = 0, x = 0onS. (11) 

The first condition has to do with the choice of the conformal factor f2, which up to now 
remains unspecified. Together they state that our investigation will be restricted to hypersurfaces 
S which are maximal with respect to both metrics g^ v and so that we can make use of the 
so-called conformal Ansatz to study the solutions of the constraint equations. Under assumption 
^ the conformal constraint equations on S reduce to 

D a D a - V) = \ XaP x aP $ with <d = rr 1 / 2 , (12) 
D a (rr 2 x Q/3 ) = o. (13) 

In order to solve in terms of expansions the constraint equations we shall make use of a slight 
variation of the so-called conformal method consisting of the following steps: 

(i) choose a negative definite metric h a p on a 3-dimensional, orientable, smooth, compact 
manifold S. In 5, choose a point i and set S = S \ i. 

(ii) Find a symmetric, ft,-tracefree tensor field ifiap on S satisfying 

D a iP a p = 0. (14) 

A standard way of finding the tensor ip a p is by means of a York splitting. That is, choose 
a symmetric, /i-tracefree tensor ip' a/3 on S and set 

2 

T/Wi = + D a vp + Dpv a - -hapDjV 1 , (15) 

where v a is some 1-form. It is well known that the latter Ansatz upon substitution on 
equation i|13|) leads to an elliptic equation for v a . 
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(iii) Finally, substitute i/) a p = d A ijj a f3 into equation (|12|l and find a positive solution d of the 
resulting Licnerowicz equation: 

(o a D a -\r\*= \^ a ^ af3 ^ 7 - (16) 

The fields h a p, Q, — d 2 and Xafi — Q 2r >Paf3 thus obtained provide a solution to the conformal 
constraint equations (|12|) and 

The above procedure is to be supplemented with asymptotic (boundary) conditions for the 
fields d and i^afi consistent with the conditions JIJ and (J5J) ■ Consider a > small enough so that 
the h— metric ball B a (i) centred on i is a strictly convex normal neighbourhood of i, and let {x a } 
be normal coordinates with origin at i based on an h— orthonormal frame {e^} at i. Consistent 
with the conditions and (|6a(l - (|Hc|l we shall require that 

|a:|i?-^l, as |a;j — »• 0. (17) 

The tensorial fields x a p, Xap, ipa/3 are related by 

Xq/3 = Vl^Xaf} = ^"0a/3 (18) 

Thus, the appropriate behaviour for the unphysical fields is given by 

Xa p = 0(l), ^ = o(^-] as |x| ^ 0. (19) 



In contrast to the analysis carried out in |15| . here we shall consider initial data sets — i.e. 
solutions to the conformal constraint equations — which are not necessarily time symmetric, that 
is, in general we shall assume that Xa/3 ^ 0. Discussions given in |2J E23 1231 wm arise as 
particular cases of our present analysis. This first article shall be concerned with the analysis of 
conformally flat initial data sets. Accordingly, we make the following assumption 

Assumption 2 (conformal flatness). It will beassumed that there is a neighbourhood, B a (i) C 
S of radius a centred on i and coordinates {x a } for which the (unphysical) conformal metric is 
of the form 

h a b = —S a b = diag{— 1,— 1,— 1). (20) 



Results regarding the nonexistence of conformally flat slices imply that the latter assumption 
eliminates all stationary initial data sets except (time asymmetric) Schwarzschildean ones from 
our considerations [201 1251 — this observation will play an important role in the interpretation 
of the results of our investigations. Moreover, based on the results of Dain [71|S] it seems that any 
discussion of strictly stationary data — i.e. data whose development is stationary but not static — 
has to consider conformal metrics which are non-smooth. For example, the conformal metric one 
obtains for Kerrian initial data form t = constant slices in Boyer-Lindquist coordinates is just 
C 2 ' a at infinity 4 . This non-smoothness of the conformal metric shall be discussed in extensis in 
paper II. 

2.1 Existence of solutions to the conformally flat constraint equations 

Under the assumptions Q and J2J) the Einstein constraint equations reduce — in Cartesian 
coordinates — to 

d a r h - 0, (21a) 

Atf = ^abV b ^\ (21b) 

4 / £ C' p ' a means that the function has p-th order derivatives f which are Holder continuous with exponent a. 
A function / is said to be Holder continuous with exponent a at a point xq if there is a constant C such that 
\f[x) — f(xo)\ < C\x — xo\ a , < a < 1 for x in a neighbourhood of xq. 
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together with the boundary conditions l|17Jl and 119(1 . We want to consider solutions to these 
equations which can be expanded in powers of \x\ — note that f(x a ) = \x\ is not a smooth function 
of the coordinates x a as d a f = — x a /\x\. In particular, we shall exclude from our discussion any 
solutions whose expansions contain fractional powers of |a;| or hi\x\ terms. These considerations 
lead naturally to the following definition, which we retake from |SJ: 

Definition (E°° spaces) . A function f 6 C°°(5) is said to be in E cc (B a (i)) if on B a (i) one can 
write f = h + \x\h with h, h e C°°(B a {i)). 

The type of solutions of the constraint equations that we want to consider are naturally 
described in terms of these spaces. The existence of the solutions to ((16(1 and ((14(1 are guaranteed 
by the following two results due to Dain & Friedrich 9 . Note that these results apply to a 
larger class of data than the one under consideration — here h ab = S ab , which is trivially smooth. 
However, these are not general enough to encompass strictly stationary initial data. 

Theorem 1 (Dain & Friedrich, 2001). Let h ab be a smooth metric on S with non-negative 
Ricci scalar r. Assume that i/} ab is smooth on S and satisfies on B a (i) 

\x\ s ip ab r b eE°°(B a (i)). (22) 

Then there exists a unique solution of equation i(_?6() which is positive, satisfies the boundary 
conditions Ji7| ) and \iy\) and in B a (i) has the form 

0=O, ^e-E°°(B a (i)), <?(*) = 1. (23) 
\x\ 

The condition \x\ s 4'a,bi) ab € E°°(B a (i)) excludes the possibility of second fundamental forms 
with linear momentum — the so-called boosted slices. In [Sj it has been shown that if one considers 
a smooth conformal metric and second fundamental forms with linear momentum, then the 
solution of the Licnerowicz equation will contain logarithmic terms at order \x\ 2 . 

Theorem 1 is complemented by the following result stating the conditions for the existence of 
a second fundamental form of the required type. It reads 

Theorem 2 (Dain &: Friedrich, 2001). Let h ab be a smooth metric on S. There exist traceless 
tensor fields ip ab € C°°(S \ {i}) satisfying equation with the properties: 

(i) in the normal coordinates {x a }, 

A 3 
?Pab = ,-|3 (3n a n b - 6 ab ) + —^(n b e cad J d n c + n a e dbc J c n d ) + 0{\x\' 2 ), (24) 
\x\ \x\ 

with n a = x a /\x\ . 
(ii) D a i; ab = on S 
(m) \x\^ ab r b eE°°(B a (i)). 

The form of the leading terms of the extrinsic curvature given in the latter theorem is, except 
for the absence of a term containing linear momentum, that of the Bowen-York Ansatz, which 
is routinely used in numerical simulations of black hole collisions — see e.g. Q]- Note that 
the theorem allows for the presence of higher order terms — which, as it shall be shown in the 
sequel, will encode the multipolar content of the second fundamental form. These terms will play 
a crucial role in paper II when discussing the structure of stationary data. In relation to this 
point, we note that the results of |21] strongly suggest that the only stationary data contained 
in the hypothesis of theorems 1 and 2 are non-time symmetric Schwarzschildean data. We shall 
return to this point later. 
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3 The initial slice near i 



Consider the Cartesian coordinates {x a } discussed in the previous section, and the orthonormal 
frame {e^} such that — x a /\x\. Under the assumptions made in the previous section one has 
that the solutions of the constraint equations l|14|l and lead to a conformal factor of the form 

n = \x\ 2 -m\x\ 3 + (D(\x\) 4 . (25) 

Let d l j kl = Qr 1 C % j k i be the rescaled Weyl tensor, where C % ^ kl is the Weyl tensor of the metric 
g^u. The expansion (|25Jl implies for the components of d l j kl an expansion in B a (i) of the form: 

m(h lk - 3<5?<M<5° /l 



where hij = h(ei, ej) — habC^e^ = —Sij. In order to discuss this kind of singular quantities, it is 
convenient to consider a certain submanifold of the bundle of frames. We present here the crucial 
points of the construction of this manifold. For full details the reader is remitted to ^S]> QB) or 

[Ej. 

3.1 The manifold C a 

In what follows we shall be using a space spinor formalism analogous to a tensorial 3+1 de- 
composition 5 . Consider the (unphysical) spacetime (.M, <? MJ/ ) obtained as the development of the 
initial data set (<S, h a p, Xap)- Let SL(S) be the set of spin dyads S — {Sa}a=o,i on S which are 
normalised with respect to the alternating form e in such a way that 

e(5A,8 B ) = £ab, £oi = 1- (27) 

The set SL(S) has a natural bundle structure where S is the base space, and its structure group 
is given by 

SL(2, C) = {t A B e GL{2, C) | e AC t A B t c D = e BD }, (28) 

acting on SL(S) by S >— > S ■ t = {&At A b\b=q.i- Now, let r = \/2eo, where eo is the future g-unit 
normal of S and 

taa> = 9(r, 6 a 6a>) = e A °€ A ?' + e A x e A , V (29) 

is its spinorial counterpart — that is, r = r l ei — o 1 aa ,t aa where a\ A , denote the Infeld-van 
der Waerden symbols and {ei}i=o,...,3 is an orthonormal frame. The spinor t aa > enables the 
introduction of space-spinors — sometimes also called SU(2) spinors, see PE1E2 : it defines a 
subbundle SU{S) of SL(S) with structure group 

SU(2, C) = {t A B e SL(2, C) | T AA d A B t A ' B , = t bb ,}, (30) 
and spatial van der Waerden symbols: 

„AB (A B)A' i ^ A'i i o mi\ 

a t =(T i A' T ' a AB= T (B a A)A>i « = 1,2,3. (31) 



The latter satisfy 



h%j — O-iAB&j , —SijCF AB a J CD = —£A{C e D)B = h-ABCD, (32) 



with hij = h{ei, ej) = —5^. The bundle SU(S) can be endowed with a su(2, C)-valued connection 
form lo a b compatible with the metric h a p and 1-form a AB , the solder form of SU(S) 6 . The solder 
form satisfies by construction 

h=h al3 dx a ®dx f) = h AB CDO- AB ®o- CD , (33) 



5 Howover, it must be said that the use of such a spinorial formalism is not essential. An equivalent discussion 
can be carried out in terms of frames. See e.g. 1161 

6 To be more precise, these structures are inherited from their analogues on 0+(S), the bundle of positively 
oriented orthonormal frames via the covering map of SU{2) onto the connected component 50(3) of the rotation 

group given by SU(2) 9 t A B ► t ! . = ° AB t A c tB D a f D e s °(3), the induced isomorphism of Lie algebras , #« 

and suitable contractions with the Infeld-van der Waerden symbols. 
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where a AB = a AB dx a — note that the o AB are not the spatial Infeld-van der Waerden symbols! 

Now, given a spinorial dyad <5 G SU (S) one can define an associated vector frame on + (<S) 
via ej = ej(S) — o ab 8at b B S B >. It is noted that different dyads 8 and 8' can give rise to the 
same frame as long as they are J7(l)-related 7 . We shall restrict our attention to dyads giving rise 
to frames {ej}j=o,--- ,3 on B a (i) such that e% is tangent to the /i-geodesics starting at i. Let H 
denote the horizontal vector field on SU(S) projecting to the aforcdiscusscd radial vectors e%. 

The fiber 7r —1 (z) C SU{S) — the fiber "over" i — can be parametrised by choosing a fixed 
dyad 8* and then letting the group SU(2, C) act on it. Let (—a, a) 9 8(p, t) e SU(S) be the 
integral curve to the vector H satisfying 8(0, t) = 8(t) e 7r _1 (i). With this notation we define the 
set 



C a = ^S(p,t)eSU(S) | |p|<o, t€SU(2,C)j (34) 

which is a smooth submanifold of SU(S) diffeomorphic to (—a, a) x SU(2,C). The vector field 
H is such that its integral curves through the fiber 7r _1 (i) project onto the geodesies through i. 
From here it follows that the projection map ir of the bundle SU(S) maps C a into B a (i). 

Let in the sequel T° = 7r _1 (i) = {p = 0} denote the fiber over i. It can be seen that 
T° w SU(2,C). On the other hand, for p £ B a (i) \ {i} it turns out that 7r _1 (p) consists of an 
orbit of U(l) for which p = \x(p)\, and another for which p = — \x(p)\, where x a (p) denote normal 
coordinates of the point p. Because of the latter, it is convenient, in order to understand better 
the structure of the manifold C a to quotient out the effect of U(\) on C a (i). The projection 
map restricted to C a C SU(S) can be factorised as C a C' a = C a /U(l) B a (i). Now, 
the set TT2 1 (B a (i)) consists of a component on which p > and another, C'~ on which 
p < 0. The projection tt2 maps both C'^~ and C'~ into the punctured disk B a (i) \ {«}. This 
last fact can be used to identify B a (i) with acquiring in the process a boundary given by 
7ri(2°) = n^ii) « S 2 . The set B a (i) = (B a (i) \ {i}) U ^{i) ~ [0, a) x S 2 is a smooth manifold 
with boundary. Summarising, we have obtained an extension of the initial (physical) manifold S 
by blowing up the point the point i into a sphere. 

For practical purposes, it is more convenient to make use of the 4-dimensional ?7(l)-bundle 



Ct =C+(J1° = {S = S(p,t)eC a | p{8) > 0\ « [0,o) x 5(7(2, C). (35) 



The construction of the manifold C a is such that a number of useful structures are inherited by C a 
from SU (S). In particular, the solder and connection forms can be pulled back to smooth 1-forms 
on C a - These shall be again denoted by <r AB and u) A B . They satisfy the structure equations: 

da AB = -uj a e A o EB - 6j b e A o ae , (36a) 
duj A B = -lu a e A lu e b + n A B , (36b) 

with 

il B = 2 r BCDEF a Acr i 

the so-called curvature form determined by the curvature spinor tabcdef given by 

TABCDEF = {^SABCE - J^rtlABCE^J ^DF + (^SABDF - -^rhABDF^JeCE, (38) 

where sabcd — S(abcd) is the tracefree part of the Ricci tensor of h a p and r its Ricci scalar. 
These satisfy the 3- dimensional Bianchi identity 

D AB sabcd = \d C dt. (39) 



7 Here we shall use the realisation 

[/(!) = {t G 1/(1) | * = ( e _° i( 
Thus, <5 and 5' are said to be [/(l)-related if 5' A = t B A 5 B with t = (t A B ) e (7(1). 
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In the case of the data considered in this paper — conformally flat, assumption 2 — it follows 
that sabcd — and r — so that Cl A B — 0. A detailed analysis of the solutions of the structure 
equations in the case of non-smooth conformal metrics will be given in paper II. 

3.2 Lifts to C a and reality conditions 

One can use the function p and the matrices t — (t A B ) G SU(2,C) to coordinatise C a . Again, let 
{x a } denote normal coordinates on B a (i) such that x a (i) = 0. The projection map ir oiC a onto 
B a (i) has the local expression 

tt : (p, t) i ► x a (p, t) = V2MX^ C cA- (40) 

This last expression can be used to "lift" scalar and spinorial fields on B a (i) to C a . In an abuse 
of notation we shall, generally, denote the lift to C a of a spinorial field pab---ef on B a (i) again 
by Hab—ef- To which manifold the field belongs will be clear from the context. 

The following two maps between spinors will play a role in the sequel: given a "primed" 
spinorial field fiA', one can associate an unprimed one via \ia : l_ * HA = t a a ha'- The other one 
is the Hermitian conjugation map fiA i— » fi \ = t a a ~p A , . These maps are extended to higher 
valence spinors in a direct way. A spinorial field VA 1 A' 1 ~-A k A' k is said to be spatial if and only if 
T A; > A iUA 1 A' 1 ---A j A' j --A k A' k = for j = 1, . . . , k. In which case its "unprimed" version is given by 

HA 1 B 1 -A k B k = T B x Al ■ ■ ■ T B k A '" v A 1 A' 1 ---A k A' k = ^(A 1 B 1 )---(A k B k ) ■ (41) 

Conversely, any spinor field with 2k unprimed indices corresponds to a spatial spinor. Thus, for 
a given spinorial field haa> one obtains an orthogonal splitting of the form 

fJ-AA' = TCA'T HAB' = -jTAA'T HCC ~ T A 'HCA, (42) 

where hca = T (c' B l jL A)B l corresponds to the spatial part of Haa'- Again, this procedure can 
be extended in a standard way to higher valence spinors. As a final remark we note that, 
given a spatial spinor field HA 1 B 1 --A k B k = H(A 1 B 1 )---(A k B k ) the frame components Hji---j k = 
o~ AlBl ■ ■ ■ (J A k kBk HA 1 B 1 ---A k B k correspond to those of a real tensor field if and only if it satisfies the 
reality condition 

l J -A 1 B 1 -A k B k = (-1) MAiBi...A fc B fc . (43) 

3.3 Decompositions in terms of irreducible spinors 

A symmetric valence 2 spinor /jab has 3 essential components. In order to make this explicit, we 
make use of the irreducible spinors 

x A b = s/2e {A °e B ], Dab = --^=e A 1 e B 1 , z AB = -^=e A °e B °, (44) 

to write 

HAB = HxXAB + tlyVAB + Hz ZAB- (45) 

Similarly, a valence 4 spinor vabcd = v (ab)(cd) nas 9 essential components and can be written 
as 

12 3 4 

VABCD = V ^ABCD + v l e ABCD + v 1 e ABCD + V &ABCD + u i e ABCD 

+ Vx(tACXBD + tBDXAc) + V y {e A CVBD + ^BDVAc) + V z (e A CZBD + i-BDZ-Ac) + V h h A BCD, (46) 

where e\bcd = € (a^ E e B F e c° € d)^' ■ notation, which shall be used again in the sequel, (•••), 
means that the indices are to be symmetrised and then i of them should be set equal to 1. In 
particular, if vabcd is associated with a symmetric tensor v a {i — ^( Q( 3) then v x = v y = v z = 0. 
If, furthermore, z^ Q/ 3 is traceless then V] x = 0. 
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3.4 Vector fields on C a 

Consequently with the discussion of the previous section one has that H — d p . Vector fields 
relative to the SU(2, C)-dependent part of the coordinates are obtained by looking at the basis 
of the (3-dimensional) Lie algebra su(2,C) given by 

1/0 i \ 1/0 -lA I ( i \ 

Ul = 2[i ) ' U2= 2{l ) ' " 3 = 2 { -i J • (4?) 

In particular, the vector is the generator of U(l). Denote by Zi, i — 1,2,3 the Killing vectors 
generated on SU(S) by ut and the action of SU(2, C). The vectors 2Tj are tangent to 1°. On 1° 
we set 

X+ = -(2 , 2 + tZi), X-=-(Z 2 -iZi), X = -2iZ 3 , (48) 

and extend these vector fields to the rest of C a by demanding them to commute with H = d p . 
For latter use we note that 

[X 1 X + ]=2X +1 [X,X_]=-2X_, [X+,X_] = -X. (49) 

The vector fields are complex conjugates of each other in the sense that for a given real- valued 
function W, X_W = X + W . More importantly, it can be seen that for p E B a (i) \ {i} the 
projections of the fields H, X± span the tangent space at p. 

We define a frame cab = C(ab) dual to the solder forms a CD and require it not to pick 
components along the fibres — i.e. along the direction of X. These requirements imply 

(a AB ,ccn) = h AB CD , c CD = c l CD d p + <£ D X+ + c^X_. (50) 

Let and a be 1-forms on C a annihilating the vector fields d T , d p and having with X± the 
non- vanishing pairings 

(a + ,X+) = (a-,X-) = (a,X) = 1. (51) 
From the properties of the solder form a AB one finds that 

C AB = X AB, c\ B = - p ZAB + c\~ B , C^ B = -VAB + C AB - (52) 

In particular, for the case under consideration — assumption 2 — , one has c^ B — 0. Furthermore, 
one has that 

a AB = ——x AB dp - 2y AB a+ - 2z AB o~ , (53) 
P 

so that the projection of (|33|l onto B a (i) renders 

h = -dp ®dp- 2p 2 {a + ® a~ + a~ (g> a + ) = -dp 2 - p 2 da 2 , (54) 
so that da 2 — 2(a + ® a" + or ® a + ) corresponds to the pull back of the standard metric on S 2 . 

3.5 The connection coefficients 

The connection coefficients are defined by contracting the connection form & A B with the frame 
cab- I n general, we write 

lCD A B = (U A B,CCD) = - p l*CD A B + lCD A B, (55) 

where, 

lABCD = 2 i e AC x BD + tBDXAc), (56) 

denotes the singular part of the connection coefficients. The regular part of the connection can 
be related to the frame coefficients cab via commutator equations which under assumption 2 are 
trivially satisfied as Jabcd = 0. We defer any further discussion of these matters to paper II. 



11 



Let / be a smooth function on B a (i)- We denote again by / its lift. The covariant derivative 
of / is then given on C a by 

D AB f = c AB {f). (57) 

Similarly, let /iab represent both a smooth spinor field on B a (i) and its lift to C a - Then the 
covariant derivative of \xab is given by 

DabUcd = cab(^cd) - 1ab E c^ed - 1ab E d^ce- (58) 
Analoguos formulae hold for higher valence spinors. 



3.6 Normal expansions 

In order to study the behaviour of diverse fields on C a near 1° in a detailed manner, we shall make 
use of a certain type of expansions which are obtained by lifting the Taylor expansions (along the 
radial direction) of the fields and then conveniently symmetrising. 

A notion which will be of much help is the following: a smooth function / on C a is said to 
have spin weight s if 

Xf = 2sf (59) 

where 2s is an integer. If furthermore, / is analytic in B a (i) then on C a it admits an expansion 
of the form 

oo 2q 

/=E^E EW^V^' (60) 



»! 

p=\s\ q&Q p m=0 



where 



Q p = |q € Z + U {0} | \s\ <q<p, q even if p is even, q odd if p is odd 1 (61) 

and fp-2q,m £ C. The functions T m J k appearing in the expansion H60[) are the functions of SU(2, C) 
onto C given by 



1/2 , \ 1/2 



SU(2,C)3t~Tj k (t)=( m ) r) t^W-t^Vv. ( 62a 



T ° (t) = l, j,k = 0,...,m, m= 1,2,3,..., (62b) 

where the string of indices with a lower index k, say, means that the indices are symmetrised and 
then k of them are set equal to 1, while the remaining ones are set equal to 0. The function / 
will be said to have axial symmetry if its expansion l|60|l is of the form 

oo 

/=E^E wv^ p - ( 63 ) 

p=\s\ q&Qp 

The following properties of the functions T m 3 k are crucial: under complex conjugation they 

transform as 

T m \ = (-iy +k T m m -\ n _ k . (64) 

Furthermore 



X+T m k } = Vj(m-i + l)T m *_! XJT m k } = - V(j + l)(m-i)T m * +1 , 
AX m fc , = (m-2j)T m fe r 



Also important for our aims is the fact that the set {\/m + IT J k } constitutes an orthonormal 
basis of the space L 2 (fi, 517(2, C)) where /x is the Haar measure on 5f/(2,C). 
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The functions T m k j are closely related to the spherical harmonics. These are defined as a 
system of orthogonal functions on S 2 . They can be extended to £> 3 as functions with zero spin 
weight. From this it follows that they can be expanded as Yi m = ^ c kjT 2k J k . Indeed, one has 

s Y nm = Hr +2 "-" l ^^T 2n "-™ _ s . (65) 

Any analytic function on SU (2, C) can be expanded in terms of functions T 1 J k . In particular, 
the product T 2ni kl ti x T 2n ^ 2 l2 can be linearised rendering 

711+712 

T 2 m \ xT 2n 2 \ = E (-l)" + " 1+ " 2 C(«i,«i-'i;«2,n 2 -^2;n,n 1 +n 2 -/ 1 -/ 2 ) 

n=q 

xC(m, ni - h; n 2 , n 2 - k 2 ;n, n x + n 2 - h - k 2 ) x T 2n n+k ^- n ^\ +h+h _ ni _ n2 , (66) 

where Qq = max{|n! — n 2 \, ni + n 2 — ki — k 2l n\ + n 2 — — l 2 } and C{1\ 1 mi] l 2l m 2 ; Z, ni) denote 
the Clebsch-Gordan coefficients of SU(2, C) 8 . 

Finally, we note that in occasions we will be in the need of dealing with fields on B a (i) which 
do not project to analytic fields on B a (i), but still have a well defined spin weight. Let g be one 
of such fields with: 

9 = E ( 67 ) 

p=0 1 

where 

n(p) 2q 

5p= E E 9 P :2 q , m T 2q m q _ s , (68) 

qr= | s | 777=0 

then we say that the coefficent g p has expansion type n(q). 

We finish by mentioning that the components of a given spinor arising in the decompositions 
(|45|l and 146fl in terms of elementary irreducible spinors have definite spin weights. More precisely, 
given the rank 2 spinor (iab, its components /iy, [i x and \x z have spin weight —1, and +1 
respectively. While for a rank 4 spinor, vabcd, the components v 2 , vu and v x have spin weight 
0; v\, v z have spin weight +1; v% and v y have spin weight —1; while vq and z/4 have spin weight 
+2 and —2 respectively. 



4 Expansions of solutions to the momentum constraint 

Using the framework described in the previous section, we proceed now to discuss the solutions 
to the conformally flat momentum constraint. The essential content of this section is not new, 
but the presentation is. Other approaches can be found in, for example, 121011 • 

We begin with some fairly standard considerations. We shall refer to R 3 endowed in standard 
(Cartesian) coordinates {x a } with the metric — 8 a b as to the Euclidean space E 3 . In 9 it was 
shown that any smooth solution in B a (i) \ {i} C E 3 of the momentum constraint can be written 
in the form 

r b = v-p + rj b + n b + r Q b + v-f , (69) 

where tpf = ©(M" 4 ), ipf = 0{\x\- 3 ), Va = ^(\x\- 3 ), and = 0(\x\- 2 ) correspond, re- 
spectively, to the parts of the second fundamental form associated with translational, rotational, 
expansion, and boost conformal Killing vectors. They all consist of I = 0, 1 spherical harmonics. 
Let n a = x a /\x\ : then if ij}f ^ the vector 

P« = J_li m / \x\ 2 ^ bc n b (2n c n a - 5 ca )dS e , (70) 

8 Some other notations used in the physics literature are: 

C(2l, 7711 ; ^2 , 7712 ; I, 777 ) = (Zl , Z2; 7711 , 7772 K, 777 ) = CQ±, 1%, l\mi,m2, 77l). 
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corresponds to the ADM linear momentum of the initial data. Similarly, if 0j fc 7^ then 

ja = _L l im f \ x \^ bc n b e cad n d dS e , (71) 

is the ADM angular momentum of the data. In the last two expressions S e denotes a sphere of 
radius e centred at i, dS e its volume element, and n a its unit normal. 

The term 0^ b in equation (|69[) is, on the other hand, associated with / > 2 spherical harmonics 
and, thus, can be interpreted as providing the higher order momenta content of the second 
fundamental form. It consists of derivatives of a certain spin-weight 2 potential 3 2 A, where 
A is an arbitrary complex C°° function in B a (i) \ {i}, and 3 denotes the standard NP "eth" 
operator. Important for our purposes is that if P a = and |x|A € E°° (B a {i)) — in which case 
tpf = 0{\x\-~ 2 )— then \x\ 8 ip ab ip ab G E°°(B a (i)) —see theorem 15 in reference 0— so that 
theorem 1 and 2 apply. 



4.1 Spinorial version of the momentum constraint 

In the space spinor formalism, a symmetric tracefree tensor Q/ 3 is represented by the totally 
symmetric spinor iPabcd = ip(ABCD) living on C a . The spinorial and tensorial versions are 
related in a standard way via 

^ABCD = OAB a CDi>U = ^AB^CD^^j^ap- (72) 

The symmetries of i/jabc d together with its tracelessness — j4£ ' AB = — imply a decomposition 
in terms of elementary spinors of the form: 

0ABCD = XO^ABCD + ^ABCD + 1>2*ABCD + ^3 e ABCD + ^ABCD- ( 73 ) 

It can be verified that the components 0o, 0i, 02, 03 an( i 04 are respectively of spin weight 
2, 1,0, — 1, —2. The components of iPabcd, being those of a spinor associated with a real spatial 
tensor, must satisfy in virtue of the reality condition (|43[) . the relations 

02 = 02, "01 = -V>3! V'O = 04" (74) 

The momentum constraint can be rewritten as 

D AB i>ABCD = 0. (75) 



4.2 The 1 — 0,1 harmonic solutions 

Arguably, the most important solutions to the Euclidean momentum constraint l|75|) are those 
having I = 0, 1 harmonics. The importance of this class lies in the fact that they are associated to 
conformal Killing vectors of the Euclidean space. Because of the spin weights of the components 
of ipABCD, we can write — cfr. the expansion Ij60|l — 

00 = O, 

01 = ^1,2,0^2 °0 + 01,2,1^2 + "01,2,2 T 2 \, 

02 = 02,O,o7o °0 + 02,2,0^2 \ + 02,2,1 T 2 \ + 02,2,2^2 \i 

03 = 03,2,0^2 ° 2 + 03,2,1^2 \ + 03,2,2J 1 2 2 2 , 

04 = 0. 

The substitution of the above expressions into equation l|75|l renders a solution which can be 
written as — cfr. equation (|69|l : 

IpABCD = tpABCD + VabCD + ^ABCD + ^PaBCD: (76) 

where the non-vanishing components of 4>abcd are gi ven by 

02^ = -^Vo- (77) 
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Those of 4>abcd are 



< = j,{P2 + iPx)T 2 ° - ^(P 2 - iP,)T 2 \, 

^ = w {P2 + tPl)T2Di " W {P2 lPl)TA 

^ = --,{P2 iPi)T 2 \ - ^P,T 2 \ + \{P 2 + lPl) T 2 % 
P P P 

While in the case of iPabcd one nas ' 

V? = \{Q2 + ^Ql)T 2 - ^fQ 3 T 2 1 -\{Q 2 - iQi)T 2 \, (78a) 

pZ pZ pZ 

^ = W iQ2 + lQl)T2 ^ - J2 Q ^ ~ ^ (02 ~ ^iWl ( 78b ) 

V> 3 Q = -\{Qi-iQi)T 2 \ - ^Q 3 T 2 \ + l(Q 2+l Q 1 )T 2 ° 2 . (78c) 

pZ pZ pZ 



And finally for ip J ABCD , 

1>{ = ^ (- Ji + iJz)T 2 ° + J 3 T 2 \ ~^(Ji + iJi)T 2 2 , (79a) 
^ = (79b) 
= p\(Ji + i-h)T 2 2 2 ~ ^i.hT 2 \ - A(_ + zJ 2 )T 2 ° 2 , (79c) 

In the above expressions A, Pi, P 2 , -P3, ^1, ^2, J3, Qi, Qi, Q3 € R. 

4.3 Solutions with higher harmonics 

The solutions to the Euclidean momentum constraint with harmonics / = 0, 1 constitute es- 
sentially what is known as the Bowen-York Ansatz [3]. This Ansatz does not exhaust all the 
possibilities as it excludes solutions with higher harmonics. Some initial data sets of containing 
these higher harmonics have been recently considered — see ^U]. The interest of these higher 
harmonic terms lies in the fact that — as will be discussed in detail in the article II — the second 
fundamental forms of stationary data, and in particular the Kerr data, generically do contain 
this kind of terms. An example of conformally flat initial data with higher order multipolcs 
can be found in ^U]. In that reference, the second fundamental form is essentially that of the 
t = constant slices of the Kerr spacetime in Boyer-Lindquist coordinates. 

Retaking the strategy used for the / = 0, 1 harmonics solutions, and because of the linearity of 
the equations we shall consider the following Ansatz for the solutions of the momentum constraint: 

2q 

V>« =^2^L n , 2q . k {p)T 2q k q ^ 2+nl (80) 

q>2 k=0 

where n — 0, . . . , 4. Substitution of the latter experssion into equation (|75[1 yields the following 
system of ordinary differential equations for the coefficients Lx-2 q ,ki L 2 . 2q ^ k and Ls-2q,kj 9 > 2, 
k : = 0, . . . , 2q: 

t . . . 1 3 2 

[q + l)L 2 -.2 q .k + -L^k - -V(? + 2)(?- l)Lo. 2g , k = 0, (81a) 

P P 
3 3 

[q + l)i 3 ;2g,fe ~ —L 1; 2q,kVl(l + ^)+ -^2:2q.k = 0, (81b) 



dLx-2q,k 


1 


dp 




dL2-2q,k 


3 


dp 






i 


dp 


~ 3P 



y/q{q + l)L 2 , 2q . k + -L 3[2qtk - -y / (q + 2)(q-l)L 4 , 2q . k = 0. (81c) 
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Thus, one sees that if one provides the functions Lo;2g,fc = ^0;2q,fe(p) and L^ 2q ,k = ^4:2<j,fc(p), 
then it is possible to integrate the above equations to obtain Li : 2 q ,k, L 2:2q ,k and L^q.k- In order 
to fulfill the reality conditions lf74l one has to require L -2 q ,k = (— ^) q+k L^q^k — then, the reality 
conditions of the other components will be automatically satisfied. In particular, one could take 
a complex function A £ C°°(B a \ {i}) and require 

2q 2q 

E E L w(p) = X+X+k, J2 E L ^AP) = X-X-A. (82) 

q>2 k=0 q>2 k=0 

The operator X + corresponds essentially to the 3 operator of Newman-Penrose — see ^H). The 
combination |;r| 4 A can be identified with the function A in section 4 of HJ. In order to satisfy the 
hypothesis of theorem 1 and 2 in section 3.1 — cfr. theorems 15 of [5] — we require that 

\x\ 4 A e E°°(B a (i)), i.e. A=-^ + -^, M , v e C°°(B a (t)). (83) 

\x\ \x\ 

Note that under the above requirements X + X + p = 0(p 2 ) and X + X + v = 0(p 2 ). Thus, one has 

oo n 2q oo n 2q 

^ = -4 E E E L^ w T 2q \_ 2 p n , V^4 = -J E E E ^ n - W T 2q k q+2 p n . (84) 

" ri=2 q=2 k=0 ^ n=2 q=2 k=0 

Under these assumptions the solutions to the system (|81all - (|81cll can be readily calculated yielding 

oo n 2q 

I'' 



oo n 2q 

^ = -a E E E ^-.n-^fc^^+i^. j = L 2 ' 3 ' ( 85 ) 



n=2 g=2 fe=0 



where 



((£4,7i-4;2g,fc - io,n-4;2g,fe)(g + 1)? + 4 J L ,r l -4;2g,fc(n- + 3) 2 ) /- — — - — 

Ll -" W = (n + 3)(2(n + 4)(n + 2)-( g + 2)( g -l)) ^ + 2){ q -l), ^ 

j- 3(L ,n-4;2<j,fc + ^4,n-4;2 9 ,fc) /7 77777 j 77 7 77 /OGU\ 

I*,n-4*,,* = 2(n + 3)2 _ g(g+1) V(g + 2)(« + !)?(? - 1), (86b) 



3,n-4;2g,fc 



((£o,7i-4;2?,fc - Li,n-i-aq,k){q + 1)? + 4L 4 ,™-4;2g,fe(« + 3) 2 ) /- — — — 

-, .ovoc ZTv — 7— -7577 777 -VW + 2 W- 1 )> ( 86c 

(n + 3)(2(n + 4)(n + 2) — (g + 2)(g— 1)) 



and 7i = 2, 3, 4, . . .. We conclude by summarising, 

Assumption 3. The spinor ipABCD associated with the second fundamental form will be assumed 
to be on C a of the form i{jabcd = fpABCD+^ABCD+^ABCD+^ABCD where ^abcd is 9 iven b V 
equation j77| ), iPabcd equations j7#a| )- |7%c| ) ; ^abcd by j79a[ )- [7$c\ ) and iPabcd by equations 
fgfl), ffip and 0fa f> - PR>($ . 

5 Expansions of solutions to the Hamiltonian constraint 

The Licnerowicz equation (|16fl is a scalar equation which can be very easily translated into the 
space spinor language rendering: 

D AB D AB ® = ^abcd^ ABCD ^ 7 . (87) 

Throughout this section we assume that iPabcd satisfies the assumption [21 and adopt the local 
parametrisation 

= X + W, (88) 
\x\ 

with W(i) = m/2, where m is the ADM mass of the initial data set [S,h a /3,Xap)- The term 
l/\x\ corresponds to the Green function of the flat Laplacian A = D ab Dab- In the case of 
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non-conformally flat 3-geometries the term l/\x\ has to be replaced by an expression of the form 
U/\x\, where the function U can be calculated recursively by means of the so-called Hadamard 
parametrix construction — this again, will be discussed at lenght in part II. Thus, the term l/\x\ 
contains information about the local geometry near infinity. On the other hand, the function W 
contains information which is global in nature — the mass in particular. 

Under the assumptions of theorem ^ the following expansion for the function W follows: 



oa p 2q 

W = m/2 + £ £ rt w va q ,kT 2q k qP v, (89) 

p=2 q=0 k=0 

with w p -2q,k G C, and satisfying the conditions w p -2 q ,k = (—~t) q+k w P :2q,2q-k so as to guarantee 
that W is real. Substitution of the above into equation l|87|l allows us to calculate an expansion 
consistent with the Hamiltonian constraint, but which not necessarily satisfies it. This is already 
quite different with what happened in the time symmetric case where any polynomial of degree 
N of the form 

N 2q 

W = m/2 + Y.H -^P^kT 2p k pP ^ (90) 

p=2 k=Q V' 

was an actual solution of the (time symmetric, conformally flat) Hamiltonian constraint 

D AB D AB d = 0. (91) 

On the other hand, in the case upon consideration in this article, the coefficients w p -2 q ,k depend, 
in general, in a nontrivial way on the parameters determining the second fundamental form. As 
an example, consider a conformally flat slice in the Schwarzschild spacetime which is not time 
symmetric. In this case ^abcd — ^abcd- ^ follows then that the Hamiltonian constraint 
implies 

AW = 4 8 mt +wr ^r + t (92) 

A solution to the latter equation can formally be written in terms of the Green function of A as: 

A 2 f 1 

W = / - rfV (QVi 

1927r7 R 3 \ x >\6\x-x>\(#* + Wy { ' 

In this specific case, it is possible to find the explicit form of the expansion — e.g. by iterations — 
up to a given order. It reads 

(94) 

For more general second fundamental forms, these similar expansions have been calculated up to 
order 0(p s ) using the scripts written in the computer algebra system Maple V. Because of their 
size, these will not be presented here. 



5.1 Some simplifying considerations 

The conformal factor $ and the symmetric tensor tp a b obtained by virtue of theorems 1 and 
2 contain some terms which are pure gauge. The freedom remaining in our setting consists 
essentially in elements of the conformal group. In particular we have at our disposition a rotation 
and a translation. We shall make duly use of them. Firstly, a rotation can be used so that the 
Cartesian coordinates {x a } are such that the ADM- angular momentum vector J° is "alligned 
with the positive z-axis" — that is J\ = J2 = 0. 

With regard to the translation freedom, this can be used to eliminate certain dipolar terms 
from the expansions of the conformal factor Consider the inversion x a = y a /\y\ mapping B a (i) 
to the asymptotic end K 3 \S!/ o (0), and consider the conformal factor <fi = \x\"d. This last function 
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satisfies, under the hypothesis of theorems 1 and 2, an equation of the form A<fi = 0(\y\ 5 ). The 
assumptions of theorems 1 and 2 — using standard arguments of potential theory — lead to 

%l M 3 \\y\ 



= 1 + ^ + ^ + ^ oil' ( 95 ) 



where d a is a constant (dipolar) vector 9 . Now, consider the translation z a = y a + c a , where c a 
is a constant vector. One has that, 

E-^-w^iw))- w-w( 1+0 (fi))- (96) 

Where it follows that 

ml / l \ 

= 1 + 2H + W [daza ~ mCaZa) + ° \W) (97) 

so that the dipolar term in the expansions can be removed if one sets c a — d a /m — "centre of 
mass". The final step now consists in a further inversion z a — x a /\x\ 2 . Thus, without loss of 
generality one can assume that 

W = m/2 + q ab x a x b + 0(\x\ 3 ), (98) 
where q a b is a "quadrupolar term" . Summarising, 

Assumption 4. The Cartesian (normal) coordinates {x a } on B a {i) are chosen so that the only 
non-vanishing component of the ADM angular momemtum is J3 — and we shall write J — J3 — 
and the expansion of i? contains no dipolar terms, that is: W — m/2 + C(|a;| 2 ). 



6 The conformal propagation equations 

In the sequel, we shall make use of the dynamical formulation of the Conformal Einstein field 
equations developed in |14II15| . This formulation of the field equations allows the implementation 
of a regular initial value problem near spacelike infinity. Besides the Levi-Civita connection V 
of the spacetime (M,7j^ u ) arising as the development of the initial data {S,hfj, v ,Xp,u)i we shall 
consider two other connections: the Levi-Civita connection V of a metric g^ v = Q 2 ^p, v in the 
same conformal class of g^', and a conformal connection (Weyl connection) 10 V, which is not 
necessarily metric. Now, given the connections V, V and V, there are 1-forms b^, f ^ such that 

V-V = S*(&) V-V = 5(/), (99a) 

V -V = 5(6- /) - 5(e _1 de), (99b) 
where for a given 1-form d^ the tensor field S(d) is given by 

S(d)^ p = 5»d p + Sfr - g w g vX dx. (100) 

The Conformal Einstein equations are essentially equations for a frame, the connection coefficients 
of V, the Ricci tensor of V, and the rescaled Weyl tensor of g^ v . 

The Conformal Einstein equations written in terms of Weyl connections allow the implemen- 
tation of a certain gauge based on the properties of certain curves called conformal geodesies 
— conformal Gauss coordinates. The conformal geodesic equations for a vacuum spacetime are 
given by 

i^V^ = -2b v x u x" + g^x^bv, (101a) 
x»\7 v b p = b v x v b^ - ^g vX bubxg^x a , (101b) 
i?Vue$ = - bu^x* - b v x v e» + g uX x»e x g^b a , (101c) 

9 It is noted that without the theorems 1 and 2, A<p = 0(\y\~ 5 ) would actually imply that the order of the 
remainder in the expansion of </> is In \x\) — see for example lemma A in the appendix of 21 . 

10 A Weyl connection is a torsionfree connection for which parallel transport preserves the causal nature. 
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where — x^{t), = b^{r), and {e^ = e^(r)}j = o,...,3 are, respectively, a curve, a 1-form and 
a frame. For q £ iS, these equations will be supplemented by the initial conditions at r = 0: 

:r(0) = g, i(0) = eo future directed and orthogonal to <S, (102a) 
K,~ 1 ttg(e j ,e k ) = r)j k , K > 0, (102b) 
= 0, 6 M = fT^O. (102c) 

The solutions to the conformal geodesic equations with these initial conditions satisfy — see e.g. 

mini— 

f(T)=eo(r), e 2 <?(e.,,e fc ) = n jk , (103) 
with a conformal factor given by 

e(r) = K - 1 rifl-r 2 4V ; 2 " =■ (104) 

The function k contains the remaining conformal freedom in this construction. The 1-form 6 M is 
such that 

(105) 



dj(T) = Gb^ = ^-2r^, K - 1 e fc (r!)^ , 



where in this last equation j = 0, . . . 3 and k = 1, 2, 3. 



6.1 The manifold -M ajK 

We introduce some useful terminology. Let CSL(2,C) = K + x SX(2,C). Further, let CSL(M) 
denote the fibre bundle with fibres CSL(2, C). Note that SL(2, C) can be regarded in a natural 
way as a subgroup of CSL(2,C), thus CSL(A4) may be identified with the set of spin frames 
{k6a}a=o,i, where {5 A } A =o,i e SL(M), k Now, C a C SU(S) C 5L(«S), from where it 

follows that C a is also a submanifold of CSL(M). On the bundle CSL(Ai) we shall be considering 
its solder form a AA , its ^-connection form u> B , and the connection form associated to the V. 
We write 

r cc , A B = (Lo A B ,ccc), r cc , A B = (u A B ,ccc), (106) 

where the the connection coefficients r cc , A B satisfy 

FcC'AB — ^CC'BA, ^CC'AB = ^CC'AB + ^CaJbC- (107) 

Because of the divergent nature of the rescaled Weyl tensor dij k i = 0(|a;| -3 ), and accordingly 
of its spinorial counterpart 4>abcd = C(p -3 ) — cfr. equation H26JI — , we shall be in the need of 
considering rescalings, 5 t— * k 1 / 2 ^, of the spinor dyads over B a (i), in order to obtain quantities 
which are regular at p = 0. A choice of the form n — pn! ', where k' is a smooth function such that 
k'(i) = 1, would do the trick. The restriction of this rescaling defines a diffeomorphism of C a to 
C a .K = K l l 2 C a . This diffeomorphism can be used to carry to C a . K the coordinates p,t A B and the 
vector fields d p , X, X + and X- defined on C a - 

The conformal factor given by equation (|104J) can be used to define the following sets: 



M a , K = \(r,q) | geC a , K , -44 <r< 441' ( 108a ) 

l=Ur,q)eM a | p{q) = 0, \r\<l\, (108b) 

r k = Ur,q)£M a | p{q) = 0, r = ±l|, (108c) 

S± = Ur,q)eMt | j6C+,,r = ±^J. (108d) 
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Let a AA denote the solder form of M.^ K . We use it to define the frame fields caa' via the 
conditions 

(<r AA ',c BB ,) = e B A e B A ', (109a) 
caa' = c AA 'dr + c AA ,d p + c\ A ,X + + c AA ,X_. (109b) 

These fields can be split in the form 

CAA' = 7}TAA'T CC CCC ~ T B A 'CAB, (HO) 

with 

t AA 'c AA > = V2d T , (111a) 

cab = r (A B c B)B , = c%d T + c\ B d p + c\ B X + + c AB X_, (111b) 

and c° AB = at r = 0. Again, the functions p, t B , T m J k are extended off C a , K by requiring that 
they remain constant, for fixed q along the curves r i— » (r, q) E K . 

The spatial part of the 1-form dj given in equation (|105|l has a spinorial counterpart <1ab = 

B' 

t< a d B } B i . The latter can be calculated in terms of the functions U and W . Namely, 

, 9 ( Ux A b- P DabU-P 2 DabW \ 

dAB = 2 <>\ WTWf )■ (112) 

6.2 The propagation equations 

Using the conformal geodesic gauge and the 2-spinor decomposition, it can be shown that the 
extended conformal field equations given in |14l I15| under the conformal Gauss coordinates de- 
scribed in the preceeding sections imply propagation equations for: 

(i) The components of the frame c AB , fi = 0, 1, ±. 

(ii) The space-spinor connection coefficients of V, Tabcd — i~ B B Tab'CD- In m] it was shown 
that in the gauge that we are using 

r AA Taa'bc = 0, T AA 'f AA ' = 0. (113) 
The latter imply the decomposition 

Fabcd = (£,abcd - X(ab)cd) — t^abIcd^ (114) 

where the spinor xabcd agrees on S with the second fundamental form of the initial 
hypersurface. 

(iii) The spinor Qaa'bb 1 associated with the tensor 

A jk = ~R( jk ) - ^V U RaVjk + ■ (115) 

where Rjk is the Ricci tensor of the Weyl connection V. In the above expression all indices 
range 0, . . . , 3. The spinor Qaa'bb' can be decomposed as 

€>AA'BB' = ®AA'BB' + ^CABfA'B' + ^AB^A'B' + ^A'B'^AB- (116) 

Now, in our gauge 

e BB 'e A A'BB' =o- (ii7) 

The propagation equations do not contain &aa' bb', but its space spinor counterpart Qac b d - 
t c a t d b Qaa'bb'- Moreover, @abc° = 0' so we can wr ite 

&ABCD — <d(AB)CD — -j^AB^>G G CD-, (H8) 
with @(AB)CD = ®(AB){CD)- 
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(iv) The components of the Weyl spinor <$>abcd — 4>(abcd), where 

dijkl — vf A 0~f B &k C a ?° ( 4>ABCD^A'B'^C'D' + <t> A' B'C D' e AB<^C D ) ■ (H9) 



The propagation equations group naturally in two sets: the propagation equations for what will 
be known as the u-quantities, v = (c^ B , £abcd, /ab, X{ab)cd, @(ab)cd, ©g G cd)> M = 0,1,± 

9rC AB = -X(AB) EFc EF-fAB, (120a) 
drC a AB = ~X (AB ) EFc EF, (120b) 
d T £,ABCD = -X(AB) EF £ EF CD + —^(^ACX(BD)EF + ^BDX(AC)EF)f EF 

-^X(AB)(C E fD)E ~ \{£AC®F F BD + e BD® F^ Ac) - iQ^ABCD, (120c) 

drfAB = -X {A B) EF fEF + ^6/ AB , (120d) 

drX(AB)CD = -X(AB) EF XEFCD ~ @(CD)AB + &VABCD, (120e) 
EF, 



d T Q(AB)CD = -X(AB) EF ®(CD)EF ~ d T Qr]ABCD + iV2d E \ a Hb)cde, (120f) 

I 

((AB) 



dr& G G AB = -X {AB ) EF ®G G EF + V2d EF n AB EF, (120g) 



where 



VABCD = ^ {(pABCD + 0~abCd)> VABCD = -^{4>ABCD ~ ^ABCLl)' ( 121 ) 

denote, respectively, the electric and magnetic parts of 4>abcd- The quantities O, d T Q corre- 
spond to the conformal factor given by equation l|104|) awhile cLab is given by equation (|112fl . 
These quantities are available a priori from a knowledge of the solutions to the constraint equa- 
tions. Thus, the equations I|120afl - |l20g| ) are essentially ordinary differential equations for the 
components of the vector v. 

The second set of equations is, arguably, the most important part of the propagation equa- 
tions and corresponds to the evolution equations for the spinor 4>abcd derived from the Bianchi 
identities Bianchi propagation equations: 



(122a) 
(122b) 



(V2- 


2c° 01 )d T (f>o 4 


- 2co O (9 T 0i - 


- 2c^ 1 9 Q ,</'o 


+ 2cQ d a 4>i 




= (2roon — 


8rioio)</>o ■ 


f (4roooi - 


h 8riooo)0i — Q^oooo4>2, 


V2d T C 


bi - c\ x d T 4>Q 


-1- CQ d T 4>2 


- Ci X d a <t>o 


+ Co <9 a </>2 




= -(4r mo 


+ fn)(f>o + 


(2roon + 


4rnoo - 2/ O i)0i + 3/oo02 - 2r ooo^3, 


V2d T c 






- c xl d a (t>i 






= -rim</io 


— 2(rnoi - 


Yhi)<t>x + 


3(rooil + rnoo)</>2 — 2(roooi — jbo)03 


V2d T C 


63 - C xl d T (j>2 






+ CQ Q d a (t>4 




= — 2rim^ 


1 - 3/11^2 


+ (2rnoo - 


j- 4r on + 2/ O i)03 - (4r ooi - /oo)</>4, 


(\/2 + 


2c° 01 )d T fa - 


2c a n d T fo H 


- 2cQ X d a ^i 


- 2c xl d a 4>3 




= -6r m i(/) 


2 + (4rmo 


+ 8rom)( 


h + (2rnoo - 8r ioi)</>4, 



(122d) 

(122e) 
where 

^=^ABCD) V .7=0,..., 4. (123) 

The subindex j in (ABCD)j indicates that after symmetrisation, j indices are to be set equal to 
1. To the equations (|122a(l - (|122e(l we add yet another set of three equations, also implied by the 
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Bianchi identities which we refer to as the Bianchi constraint equations 

= -( 2r (01)11 - 4Tiiio)0o + (2r 00 ii - 4T( i)oi - 4r U oo)0i + 6r (ol ) OO 2 - 2r OO oo03, (124a) 
c^drcj)! - 2c o ) 1 <9 r 2 + CQ d T cj)3 + c xl d a (j)i - 2co 1 9 Q 2 + c^ d a (j) 3 

= r m i0o — (4r( i)ii — 2rnoi)</>i + 3(r on — rn o)^ — (2r ooi — 4r( O1 ) O o)03 

-r ooo^4, (124b) 
c?i9 T 02 - 2cq 1 <9 7 -0 3 + CQ d T cj)4 + c xl d a (f>2 ~ 2cq 1 9 q (/)3 + Co O (9 Q 4 

= 2r m i(/)i - er (oi)u^2 + (4r 00 n + 4r (01)01 - 2r UO o)03 - (4r 00 oi - 2r (ol ) OO )04. (i24c) 

6.3 The initial data for the conformal propagation equations 

Assume for the moment that one has a solution (n,Xap) of the constraint equations (|10a|l and 
IjlObjl . How do we calculate the initial data for the propagation equations (|12t)a|l - fl20g| ) and 
(|122a(l - (ll22e(l ? In order to do this, one has to make use of the so-called conformal constraint field 
equations — see e.g. |14) . 

Let O be the conformal factor induced by the F-gauge. We consider the further conformal 
rescaling i— > k _1 where k is a function on S such that k = k' p with n'(i) = 1. This conformal 
rescaling induces a rescaling in the spinor frame of the form 6 i— > k 1 / 2 ^. Using the procedure 
described in |14l ITK] one arrives to 

2 2 2 

Oabcd = -^D [AB D CD) Q + ^ (r + X efghX EFGH ) h ABC D + -^e A BXc D EF Def^(I^) 

<Pabcd = -^D^ AB D CD )Q. + —sabcd + -^Xef(abX C d) + ^~^ dE \aXbcd)e, (125b) 
c AB = 0, c\ B = kx ab , (125c) 

c ab = K + ^as) > c ab = k {~y AB + ^b) ' ( 125d ) 

£,ABCD = V2 (^KJABCD ~ ~(e A CK>BD + ^BD^Ac)^ , (125e) 

X(AB)CD = Kn 2 ip ABCD , (125f) 
fAB = KAB, (125g) 

with 



KAB = ( X AB d p + y~ Z AB + C.\Bj X + + \y~y AB + ~ C ABj X - J K - ( 12 6) 

In the above expressions it has been assumed that the initial hypersurface is maximal but not 
necessarily conformally flat. 

7 Transport equations at X 

The system of equations Q12UaJ) -( |r2Tig| | i (|122al) - (|122e[) . an d (|124a|) - (|124c|) allow us to introduce a 
special kind of (asymptotic) expansions in the region of spacetime near null and spatial infinity. 
Writing as before v = (c^g, £abcd, Sab, X ( Ab)cd , @(AB)Cd, ®g G cd) > /" = M,± and <P = 
(0o i 0i , 02, <fe, 04 ), then the equations H120a|l - fl20g| | can be concisely written as 

d T v = Kv + Q(v, v) + L<f>, (127) 

where K , Q denote respectively linear and quadratic functions with constant coefficients, and L 
denotes a linear function with coefficients depending on the coordinates through the functions 0, 
<9 T 0, cIab and such that L\x — 0. For the Bianchi propagation equations one can write 

V2Ed T ^ + A AB c%d^ = B(T ABCD )<t>, (128) 
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where E denotes the (5 x 5) unit matrix, A AB c AB are (5 x 5) matrices depending on the coor- 
dinates, and B(T abcd) is a linear (5 x 5) matrix valued function with constant entries of the 
connection coefficients Tabcd- On similar lines, the Bianchi constraint equations (|124a|) - (ll24c|) 
can be written as 

F AB c%d^ = H(T ABC d), (129) 

where now F AB c AB denote (3 x 5) matrices, and H (Tabcd) is a (3 x 5) matrix valued function 
of the connection with constant entries. 

The system (|127|) - (|128|l can be formally evaluated on 1 rendering an interior system of equa- 
tions which, due to L\x = is decoupled, and can be easily solved. The solutions (v^°\ 4>^°') 
to this system are universal — i.e. independent of the initial data — and can be regarded as the 
leading terms in Taylor-like expansions of the form 

"-E/^V, *~£VV. (130) 
p=0 p ' P =o 1 

The higher order coefficients in these expansions can be determined via a recursive procedure. 
To do so, we differentiate the systems (I127|l . <|128|) and (I129|l p times and then evaluate at X. The 
resulting equations are of the form 

d T v {p) = Kv {p) + Q{v (0) ,v (p) ) + Q(v {p) , v {0) ) 
p-i 

+ E (Qiv^iV^-ri) + L (J V P ~ J) ) +7>V (0) - (131) 

3=1 

From the transport Bianchi propagation equations H122af> - I|122c[l one § e ^ s 

(V2E + A AB (c%)^) d T ^ + A AB (c c AB )^d c ^ = B{Tf BCD )^ 

+ E ( P j ) {bF&cd)*™ - A AB (c%)^d^) , (132) 

where C = ±. Similarly, from the Bianchi constraint equations (|124afl - <|124c[) one obtains 
F AB {c AB )^d T ^ + F AB (c c AB )^d c 4> (p) = H(T% CD )^ 

+ E ( J ) {h{T% cd )^ - F AB {c%)^d^) . (133) 

We will refer to these equations as to the transport equations of order p. An interesting feature 
of the transport equations is that their principal part is universal — i.e. independent of the order 
p. For p > 1 they are linear differential equations for the unknowns of order p. Furthermore, note 
that the subsystem (|131fl consists only of ordinary differential equations. Now, the expressions 
(|125a|) -( |T25g| l for the initial data of the conformal propagation equations can be used to determine 
the expansion types of the quantities ti^' and <f>^ at 1°. For inital data satisfying assumptions 
1-4, it is found that they have expansion type p — 1 and p respectively. The transport equations 
(|13ip and (|132|) can be used in turn to show that the expansion type is preserved during the 
evolution. 

The transport equations l|131|) . (|132fl and l|133|) arc decoupled in the following sense: a knowl- 
edge of vd\ 4>^\ j — 0, . . . ,p — 1 together with the initial data v^\xo allows us to solve the 
subsystem (|131|) to obtain the quantities »W. With cf>( l \ k — 0, . . . ,p and I — 0, . . . ,p — 1 
and the initial data (j)^ at hand one could, in principle, solve the equations l|132|) to get <j)( p \ 
A major complication in the aforediscussed procedure is that the coefficient accompaigning the 
r-derivative in equation (|132[) is such that 

(V2E + A AB (c° AB )W) = V2diag(l-r,l,l,l,l + r), (134) 
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that is, it looses rank at the sets I — in other words, the system degenerates. Intuitively, one 
would expect this degeneracy of the system to leave some sort of imprint on the behaviour of 
its solutions on the critical sets, X^. From a partial differential equations' point of view, the 
degeneracy arises from the transversal intersection of a total characteristic of the Conformal 
Einstein field equations, the cylinder X, with standard characteristics — null infinity, ^ . 

In order to obtain a detailed analysis of the behaviour of the solutions to the propagation 
Bianchi transport equations at T^- one has to make use of the transport equations implied by the 
Bianchi constraint equations 1|133|) . 

The equations H132fl - I|133|) can be, consistently with the spin weight of its diverse terms, 
decomposed in terms of the functions k - In particular, one has 

= Y a iWl,k T 2 q k q -2+ 3 - ( 135 ) 

q=\2-j\ k=0 

with Oj,p;2g,fc = dj,p;2q,k(T). The latter strategy "reduces" the problem to the analysis of a system 
of ordinary differential equations in r. Note however, that the left hand sides of the equations 
(If 3211 and (|f 33fl involve products of the components of the vectors and <f>^ which need to 
be linearised, that is, one needs to write them as a linear combination of T- k 's according to 
formula (ffiH)) . Now, proceeding sector by sector, one can use the transport constraint equations to 
further reduce the problem to the analysis of a system of two ordinary differential equations for 
the coefficients ao^ p -2q,k and a^^q.k of 0q and ■ The remaining coefficients can be calculated 
algebraically from these two. The equations can be written in the form 

(I - r 2 K + (4 + 2(p - l)r)a' + (q+p)(q-p+ l)a = b 0iP;q ( T ), (1 36a) 

(i - r 2 K + ( - 4 + 2(p - i)tK + {q + p){q -p+ 1)04 = &Wt), (136b) 

where we have written clq and 04 instead of ao,p;2g,fc and d4 jP ;2g,fe, and where the terms bo iP;q and 
&4,p;g can be derived from the right hand sides of Ijl32|l and 1|133(1 . Their solution can be written 

as 

(137) 



\a4,p;2q,k(T)J [ H \ a 4,p;2<},fc (U)/ J Q 



where X p ^ q (r) denotes the fundamental matrix of the system (|136a|l - and the vector B p ^{t) is to 
be derived from bo, P -, q and &4,p ;g . The matrix X P) q(r) has been given explicitly in terms of Jacobi 
polynomials. More importantly, 

detX p , q (r) =/(r)(I-T 2 f- 2 (138) 

with f(r) a second order polynomial such that /(±1) 7^ 0. This means that the integrand 
in equation (|I37|I will have poles at r = ±f unless -B Pj(J (t) is very special. The main hurdle 
to be overcome in this procedure is that the explicit form of B p _ q {r) becomes more and more 
complicated as the order p of the expansions increases. Thus, in order to gain some understanding 
of its structure, one has to make use of computer algebra methods and calculate explicitly the 
solutions up to a given order. 

7.1 A regularity condition 

A detailed analysis of the expansion types of the diverse quantities involved in equations (11321) 
and (|133|) reveals that B PiP (t) — 0. Moreover, 

aO,p;2p,fc(T) = (1 - T)P +2 (1 + T)P- 2 (a ,p;2p,fc(0) + {P+1) ip P+2) Kp;2p,fc(0) - a 4 ,p;2p,fc (0)]/+ ) , (139a) 
a 4 ,p;2p,fc(T) = (1 + rf+ 2 (l - T)P- 2 (a 4 ,p;2p,fe(0) + (P+1 » P+2) Kp;2p,fc(0) - a ,p;2p,fc (0)] I- ) , (139b) 

where 

J ± =r 7 *^ _=Aln(l-r)+ Ap±2 ----- Al 



(l±r') p - 1 ( 1 T-r') p+3 (1-t)p ±2 1-r 

-B. h(l + r)+ ^ 2 + • • ■ + ^ + C, (140) 
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where the A's, -B's and C are some constants. Hence, the coefficients ao,p;2p,/c (t) and a^ yP -2 P ,k (j) 
will not be smooth at r = ±1 unless 

ao,p;2 P ,fc(0) = a 4iP;2 p,fc(0). (141) 

We shall refer to the latter as to a regularity condition. Although the original derivation of this 
result in |15| was given with time symmetric initial data in mind, the statement holds true even 
if the data is not time symmetric 11 . 

Given the aforediscussed regularity condition, the question arises whether it is possible to 
rewrite it in terms of the freely specifiable data available when solving the constraint equations 
using the conformal method. This is the content of the following theorem. 

Theorem 3. For conformally flat data such that \x\ 8, ipaf34' a S E ca (B a (i)), and independently 
of the choice of the function k appearing in the conformal factor {10$ , the following conditions 
are equivalent: 

(i) </>o, P ;2p,fc(0) = </>4,p;2 P ,fc(0), p = 2, 3, . . k = 0, . . . , 2p; 

(U) Lo,p-4:2p,fc + Li,p-4;2p,k = 0, p = 2, 3, . . k = 0, . . . , 2p. 

Proof. The result follows from looking at formula l|125b(l and counting arguments regarding 
the expansion types of the quantities involved. 

Remark. Note that the regularity condition l|141|) is trivially satisfied for time symmetric, 
conformally flat data. 

The analysis described in the present article will be concerned with initial data rendering 
solutions to the transport equations which are as smooth as possible. Accordingly, we make the 
further 

Assumption 5. The second fundamental form satisfies the condition (ii) in theorem 3. 

8 Solutions to the transport equations and obstructions to 
the smoothness of null infinity 

Prom the discussion in the previous section it follows that the origin of the non-smoothness of 
the solutions of the transport equations lies in Bianchi subsystem <|132ll . In this section we shall 
describe in a qualitative fashion the structure of the solutions to this subsystem for initial data 
satisfying the assumptions 1-5. In order to ease our calculations we choose 

k = p, (142) 

any other choice with k = pn' , where k! is smooth and such that k!{i) = 1 would not alter the 
essence of the results. 

All the results here presented have been obtained by means of scripts written for the computer 
algebra system Maple V. In what follows, we shall be systematically using the decomposition 

4 P) = E Y. a iMT 2q \_ 2+j . (143) 

g=\2-j\ fc=0 

of the components of the Weyl spinor 4>abcd- We shall use Q(r) to denote a generic polynomial 
in r, while Pfe(r) will denote a generic polynomial of degree k in r such that ^(±1) 7^ 0. It will 
be understood that <2(r) and 'Pfe(r) appearing in different equations are in principle unrelated. 

11 In 1151 Friedrich, in a tour de force, has shown that the condition <jf>o,p;2p,fc(0) = </>4,p;2p,fc(0) for p = 2,3,..., 
k = 0, . . . , 2p is equivalent to the vanishing of the Bach tensor and its symmetrised, tracefree derivatives to all 
orders at i — cfr. also with footnote 1 on the present article. 
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8.1 Lower order solutions 

The first observation that has to be made is that if the assumptions 1-5 are satisfied then the 
lower order solutions are fully regular. Indeed, 

Theorem 4. Under assumption 1-5, the solutions of the transport equations for p = 0, 1,2,3,4 
have polynomial dependence in r. Thus, they extend smoothly to the sets I . 

It should be pointed out that assumption 5 plays a crucial role here. The original analysis in 
[TS] suggests that if the regularity condition involving the Cotton- York tensor does not hold, then 
logarithmic divergences in the solutions to the transport equations can appear at order p = 2. 



8.2 Solutions at order p = 5 



As it happens in the time symmetric case, the first logarithmic divergences in the solutions to the 
transport equations appear at order p = 5. More precisely, the solutions to the p = 5 v-transport 
equations are polynomial in r. On the other hand, the solutions of the p — 5 Bianchi transport 
equations are such that 

a,,5;4,fc(r) = T 5;4 , fe (a -T) 7 -^,-(r)ln(l - r) + (1 + r) 3+J > 4 -i(r) ln(l + r)) + Q(r), (144) 
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The coefficients for the remaining sectors: a,- ;5 ,o,fc(T), aj ; 5,2,fe(T), aj;5,6,fc(j)i a J; 5 i8 , fc (r) and a j;5 , w ,k{T), 
j = 0, . . . , 4, are polynomial in r. We refer to the coefficients T^-^k as to the quadrupolar ob- 
structions at order 5. Note that although the initial data set is time asymmetric, the obstructions 
obtained at this order are time symmetric. 

The quadrupolar sector at order p = 5 is the only one rendering logarithmic divergences. More 
precisely 

Lemma 1. For initial data satisfying the assumptions 1-5, the coefficients a^^2q,k, j = 0, . ..4, 
q = 0, . . . , 5, q 7^ 2, k = 0, . . . , 2q are polynomials in r. 



8.3 Expansions at order p = 6 

Some experimentation reveals that if one uses the solutions <|144|) to the p — 5 Bianchi transport 
equations then ln(l ± r) terms will be present in the vector v^®'. Moreover, <f)^ would also 
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contain, besides the ln(l ± r) terms, expressions involving ln 2 (l + r), ln(l + t) ln(l — r) and 
ln 2 (l — t). The discussion of these non- regular solutions falls beyond the scope of the present 
investigation. Hence, in the following calculation we assume that 

T 5; 4,fc = 0, fc = 0,...4. (150) 

The latter condition could be used to solve for, say, W2-4,k- Now, at first glance it is not fully clear 
whether it is possible to have some initial data for which the obstructions vanish. As already 
mentioned in section 5, the coefficients w p - t 2 q ,k in the expansions of the function W are of global 
nature and may depend in a non-trivial way on the freely specifiable data. It may be the case 
that in order to have initial data for which the obstruction vanishes further restrictions on the 
free data are necessary. Indeed, explict calculations render the following 

Lemma 2. Consider initial data satisfying the assumptions 1-5. Let T^-^k = 0, k = 0, . . . , 4. 
Then the p = 6 Bianchi constraint transport equations are satisfied if and only if 



1023 r— / T \ 224 r-/ T T \ 3099 r- 2r „ ,,„,, 

-j^j-VOTO^o-l^fc - L 4 -i ;4 ,fc ) + — V6(^ ,0;4 : fc - i4,0;4,fc ) ~ ^-VBm L -2,4,k = 0, (151) 

for j = 0, . . . , 4. The latter implies that 

2089 



m»2;4,2 = - -^- j2 . ( 152 ) 



while the other W2 ; 4,fc 's vanish. 



If Y5;4,fc = and the previous theorem is satisfied, then the solution to the Bianchi transport 
equations are of the form: 

ai,6i4,fc(r) = T+ 4jfc (l - T?-*Vj(T) ln(l - r) + Tg ;4fe (l + t) 4 +^ 4 _,(t) ln(l + r) + Q(r), (153) 

with k = 0, . . . , 6. 

These obstructions do not have a counterpart when considering expansions of time symmetric 
initial data. We call the coefficients T^! 4 k the quadrupolar obstructions at order p = 6. Further- 
more, the obstructions are asymmetric with regard to the two disconnected parts of null infinity. 
That is, 

K,4,k = & T- 4;fe = 0, (154) 

generically. 

In order to analyse with more detail the structure of the obstructions, we shall further assume 
that the initial data is axially symmetric. In this case, only the obstructions with k — 2 arc in 
principle non-vanishing. They are given by 
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27 



The asymmetry in the obstructions can be observed in much simpler situations. Namely, 
consider initial data sets for which 

Qi = io,-l;4,2 = ^4,-l;4,2 = ^0,0,4,2 = -^4,0;4,2 = Lo, 1,4,2 = i4,l;4,2 = 0. (156) 

In that case the obstructions read 

+ _ 2198208 2 20817 2 _ 2198208 - 20817 2 

T ^- 2 ~ ~~6965~ J ~ ~\T AJ ' T ^ ~ ~696^ J " ~1T AJ ' (157) 

which clearly cannot be satisfied simultaneously. 

Now, assume one has an initial data such that Y^ 4 2 = T^ 4 2 = 0. One can use these two 
conditions to solve for Lo,-i;4,2j ^4,— 1;4,2- The latter implies that 

14 / \ ^19096 

£o,-2;4,2 = 23 ( ^0,0:4,2 — ^4,0:4,2 J — V 6 J . (158) 

The regularity condition (ii) of theorem 3 together with the reality condition l|74|) imply that the 
coefficient £o,-2:4,2 is a pure imaginary number. Similarly, because of the reality conditions, the 
first term on the right hand side is also pure imaginary, while the second is clearly a real number. 
This implies that 

J = 0, w 2 -a,2 = 0. (159) 

This result holds true also for initial data without axial symmetry. More precisely, on has the 
following 



Theorem 5. Initial data sets satisfying assumptions 1-5, if T5.4 j. — 0, k = 0, ...4 with the 

"6;. 



consistency condition of lemma\^ fulfilled, and moreover, with T^ 4 k , k = 0, . ..4 are such that 



J = 0, w 2 ;4,fe = 0, with k = 0, . . . , 4. (160) 



Following a pattern already observed in the case of time symmetric initial data, the sector 
o,,6,8,fc( r ) = will also render obstructions — which, in this case, are of octupolar nature and 
time symmetric. More precisely 

a jA;6 ,k{T) = T 6;6jfc ((1 - Tf- j V j+ i(r) ln(l - r) + (1 + r^'TV^r) ln(l + t) \ + Q(r). (161) 

In particular, one has for example that 
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V30£o, -1,6,3 + -gTj"V 301 Lo,0;6,3 — ^4,0;6,3 I + V 30 I Lo,l;6,3 — ^4,1;6,3 



(162) 



The other obstructions have a similar structure. These can be used to solve for the octupolar 
quantities u>3 ; 6,fc, k = 0, ...,6. In a similar way as to what happened when setting T^-^k = 0, 
it is to be expected that consistency conditions will arise in the expansions at order p — 7. In 
complete affinity with lemma one has the following 

Lemma 3. For initial data complying with assumptions 1-5 and such that Ts-^k = 0, k = 
0, . . . ,4, if the consistency conditions of lemma\^ are satisfied, then the coefficients a>jfi-2q,k(T) 
with j = 0, . . . , 4, q — 0, ... ,6, q =/=,2,3, k = 0, . . . , 2q are polynomial in r. 
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8.4 Expansions at order p = 7 and beyond 

The analogue of lemma [21 for the expansions at order p = 7 is the following 

Lemma 4. Consider initial data conforming with assumptions 1-5. If, moreover, T5;4 & a = 
T^ 4 fe2 = 0, and Tg ; g ^ 3 = /or fej = 0, . . . , 2i and the consistency condition in lemma is 
observed, then the p = 7 Bianchi constraint equations are satisfied if and only if 

Qi = Q2 = Q3 = 0; (163a) 
£o,0;4,fc = £4,0;4,fc, with k = 0, . . . , 4; (163b) 



85^L , ;6,fc ~ L ±A;6,kj + 241 £o,i;6,jfc - Li,i-e,k j - 303i ,-i;6,)fe = 0, with k = 0, . .. ,6. (163c) 

Remark. The conditions given in the previous lemma imply that 

Lo-2-am = 0, w 3j6ifc3 = 0, (164) 

with fcj = 0, . . . , 2i. 

If the consistency requirements given in the previous lemma are satisfied then it turns out 
that the sectors aj t 7- > o,k ( T )i QjJ&MiT), aj,7;io,fc 5 (T), a-,,7 ; i2,fc 6 ( T ) and ajj-u^r) have polynomial 
dependence in r, while the sectors o J j,7;i,k 2 { T )i a j,7-,6,k 3 ( T ) an d aj,7 ; 8.fc 4 ( T ) exhibit a distinctive 
structure. Namely, 

%7;4,fc(T) = T+ 4 fc (l - r) 9 -^,(r) ln(l - r) + T" 4 jfc (l + t) 5 +^ 4 _,(t) ln(l + r) + Q(r), (165) 

with k — 0, . . . , 4 and the obstructions T^: 4 k time asymmetric in the sense described in the 
previous subsection. They depend on A, L ,o ; 4,fc, £4,o ; 4,fc, £o,i;4,fc, £4,i ; 4,fc, £o,2;4,fe, £4,2 ; 4,fc- The 
conditions Y^ 4 k — T^ 4 k — can be used to solve for, say, £o,0;4.fc> ^4,0;4,fc- 
Similarly, one has 

a^eMr) = T+ 6k (l-T) 9 ^r 3+1 (r)Hl-r) + T- 6k a (166) 
with k = 0, . . . , 6 and T^i 6 3 also time asymmetric and depending on A, -Lo,0:6,fci £4.0:6,fc, Lo.i-6,k, 

^4,l;6,fc, £(),2:6,fc an d £4,2:6,fc- 

Finally, whe have 

%-,7;8,fe(r) = T 7;8 , fc ((1 - r) 9 -^, +2 (r) ln(l - r) + (1 + r) B+J P 6 -i(r) ln(l + r)) + Q(r), (167) 

for fc = 0, ... ,8 which renders a time symmetric obstruction depending on 1^3,6, fc, A, £o,0;8,fc> 
£4,0:8, fc, £o,i;8,fe, £4,i;8,fc) £o,2 ; 8,fc and £4,2:8, fc- The condition T 7;8ifc can be used to solve for 

W 3 ,6,3- 

As it is to be expected, requiring that the obstructions obtained at order p — 7 all vanish 
gives rise to consistency conditions on the expansions at order p = 8. In particular, one has the 
following 



Lemma 5. Assume that the hypothesis and conditions of lemma^are satisfied. If one has initial 
data such that T^ 4 fea = 0, T^g fea = 0, and Tt^m^ = 0, with ki = 0, . . . , 2i then a necessary 
condtion for the existence of solutions to the p — 8 constraint equations is that 

£o,-i,6,fc = £4,-1,6,*; = 0, with k = 0, . . . ,6. (168) 



From the results here presented it is possible to hint a general pattern in the structure of 
the solutions to the transport equations at an arbitrary order. Nevertheless, it must be said 
that a rigurous proof of such a pattern, given the intrincacy and length of the equations, would 
constitute a remarkable feat which cannot be attained with the current tools and understanding 
of the problem. 
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8.5 Asymptotic Schwarzschildean data 

The vanishing of the obstructions and the fulfillment of the associated compatibility conditions 
greatly restrict the type of conformally flat initial data under consideration. We summarise the 
results of our calculations in the following 

Theorem 6. An initial data set (S,h a /3,Xap) for which 





~~ 6AM 


1 7,4, fc 2 U > 




— L 7,6,k s 


= o, 









with ki = 0, . . . ,2i, and for which the consistency conditions in lemmas\^ ^] and[3| are satisfied 
is such that 

IpABCD = 3 tABCD + 

Remark. Initial data sets of the sort discussed in the previous theorem can rightfully be 
called Schwarzschildean up to order p = 3. Their essential octupolar (2 3 = 8) content vanishes. 
More generally, we shall say that an initial data set is Schwarzschildean to order p near inifinity 
if 

W=^ + 0(ff +1 ), ^ ABCD = -^e ABCD + 0(pP- 3 ). (169) 
2 pr 

That the solutions of the transport equations for the time asymmetric data discussed in section 
5 extend smoothly to the critical sets can be shown through an induction argument. Because 
of the manifest spherical symmetric of the set up, the only non-zero field quantities are those 
with spin- weight — that is, the only non-trivial sector is that given by T ° . It is not difficult 
to show that the fundamental matrix of the system of v-transport equations l(131() has entries 
that are polynomial in r — this result is also true even without assuming spherical symmetry. 
This implies that the solutions of the v-transport equations will be polynomial if the lower order 
terms are polynomial. More crucially is that under the current assumptions, the only non-trivial 
Bianchi propagation equation is 

Sr<t>2 = -7^X2</>2 ~ tyhfo, (170) 

where xi and Xh are, respectively, the e^ BCD and Habcd components of Xabcd — see section 
3. From here it also trivially follows that if the lower order terms in, say, the pth-order Bianchi 
transport equation are polynomial, then also cj)^ will be polynomial in r and, thus, will extend 
through the critical sets. This simple argumentation is possible only because of the spherical 
symmetry. This is to be contrasted with the case of, say, the solutions to the transport equations 
for Kerrian initial data, for which — as it will be discussed in paper II — all the components of 
the Weyl spinor 4>abcd are non-zero. 

To prove that the solutions of the transport equation for the time asymmetric Schwarzschildean 
data render series of the form 1)130(1 which are solutions to the propagation equations is a bit more 
sophisticated. It can be done using techniques similar to the ones used by Friedrich in section 6 
of |15|. In any case, th smoothness follows from the analysis of |18j . 

9 Conclusions 

Theorem and the structure observed in the expansions up to order p = 7 readily suggest the 
following generalisation of the conjecture given in [21]. Namely, 
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Conjecture. The time development of an asymptotically Euclidean, initial data set which is 
conformally flat in a neighbourhood B a (i) of infinity and for which \x\ s i['ai3'4> a € E 00 (B a (i)) 
admits a conformal extension to both future and past null infinity of class C k , with k a non- 
negative integer, if and only if the initial data is Schwarzschildean to order p m in B a (i), where 
p, = Pm(k) is a non-negative integer. If the development admits an extension of class C°° , then 
the initial data has to be exactly Schwarzschildean on B a (i). 

Remarks. Firstly, it is important to notice — as it follows from our discussion in section 8 — 
that the Schwarzschildean data is not necessarily time symmetric. The condition \x\ s ipa/3'4 !al3 G 

(B a (i)) is used here for convenience. Any condition which would ensure that the conformal 
factor arising from the Licnerowicz equation is of the form $ = with $ £ E°° — a conformal 
factor that is not necessarily smooth, but whose expansions in B a (i) do not contain logarithmic 
terms — would serve as well. One could think of generalisations of the conjecture where one 
requires that "9 £ E m {B a {i)) for some given non-negative integer m. If the condition is sharp, it 
is to be expected that it would generate its own type of obstructions to the smoothness of null 
infinity. Finally, we notice that — as seen from the results of our calculations — future and past 
null infinity could have different degrees of regularity, say C k+ and C k ~ . This extra structure is 
lost in the conjecture as clearly k = min{fc+, 

A first step toward proving this conjecture would be to show that the solutions to the pth 
order transport equations extend smoothly to the critical sets if and only if the initial data 
is Schwarzschildean to a certain order p, — p,(p). The calculations shown in this article suggest 
that at least p, = p — 4. This is to be contrasted with the results in [23] which suggest that 
for conformally flat, time symmetric data the corresponding value ought to be p. = p — 3. This 
result which should provide sharp conditions could then be used in turn to prove the necessary 
existence results required for the conjecture. 
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